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Abstract 

In view of the recent interest in the decays of mesons into a pair of light 
leptons, a computation of the QED radiative corrections to the decay of 7r° 
into an electron-positron pair is presented here. The results indicate that the 
peak value of the differential decay rate is reduced by about 50%, because of 
soft-photon radiation. The number of e + e~ pairs having invariant masses in an 
energy bin of 0.5 MeV centered around is found to be about 20% smaller 
than the one it would be if radiative corrections were neglected. 







The decay of 7r° into a pair of electrons via two virtual photons (see Fig.^J) 
was studied long ago |jj . The authors of Ref . calculated the absorptive part of 



that decay amplitude, in order to derive a unitarity bound for the branching ratio 



using as input the experimentally measured value for r(7r° — > 77). They found that 
B > 4.7 x 10~ 8 , an exact result. Similar exact unitarity bounds can be derived for 
the decays of heavier mesons into lepton-antilepton (Z + /~ ) pairs. 

To complete the above calculation though, one would also need the dispersive 
part of the amplitude. The knowledge of the exact branching ratio is essential, 
because it could provide us with a test of new physics 0. To be more precise, 
the decay of a meson to a l + l~ pair could be influenced by the existence of point- 
like effective 4- Fermi interactions between the constituent quarks of the meson and 
the final two leptons. Such interactions exist in technicolor models, in models of 
composite quarks and leptons, etc. In technicolor models, for example, extended- 
technicolor forces couple not only technifermions to ordinary fermions, but also 
ordinary fermions to themselves. 

The dispersive parts of the two photon contribution to the mesonic decay 
amplitudes are difficult to calculate, however, since in that case the two intermediate 
photons are off-shell; there exists presently no exact model describing the non- 
perturbative meson dynamics needed to make this calculation. In addition, the 
problem is even more complex in the case of the pion. This is because this meson 
plays also the role of a Goldstone-boson, arising from the spontaneous breakdown 
of the SU(2)l x SU(2)r chiral symmetry of QCD, and the assumption that it is 
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Figure 1: The decay tt — > e + e . The 4-momenta of the various particles are given 
inside parentheses. 

just a bound state of two quarks must be viewed with caution. 

There have been several attempts to estimate the dispersive part of the decay 
amplitude, using a vector dominance model ||, leading to B = 1 x 1CP 7 , a nucleon 
loop model leading to B = 1.4 x 10~ 7 , an extreme-non-relativistic model [§, ||], 
leading to B = 6.2 x 10~ 8 , a relativistic bag model ||, leading to B 1 x 10 -7 , and 
a chiral Lagrangian approach Q, leading to B = 7± 1 x 10~ 8 . Early experimental 
results [|| gave B = 1.8 ± 0.6 x 10 -7 , which was somewhat large, and caused some 
controversy over the theoretical predictions ||. Nevertheless, the latest experiments 
in Fermilab |l0| and Brookhaven [i~l| ] , which give B = 6.9 ± 2.8 x 10" 8 and B = 
6.0 ± 1.8 x 10 -8 respectively, indicate that the dispersive part of the amplitude is 
indeed smaller that the absorptive part. 

Another issue relevant to the decay 7r° — > e + e~ is the one related to final-state 
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QED radiative corrections. The fact that m n 3> m e makes the two electrons in the 
final state quite energetic, and radiative effects can therefore be quite important. 
Similar instances, where QED radiative corrections are large, have already been 
observed in different mesonic decays into leptons |jL^] . Radiation of soft photons 
can alter substantially the decay profile of the pion, reducing considerably the peak 
differential decay rate of this process, and obviously affecting the branching ratio B. 
This paper is focused on these corrections, and their treatment is done independently 
of the method used to calculate the non-radiatively-corrected decay rate. 

An early calculation of QED radiative corrections to the decay of a meson (a 
kaon) into two leptons via two virtual photons, to first order in the fine structure 



constant a, appeared in Ref.|l3|]. Later, a similar calculation appeared [14] for 
7T° — > e + e~ , the decay studied here. That calculation was also limited to first order 
in a. In such a process, however, e + e~-invariant-mass cuts must be very strict, as 
much as experimental resolution allows, in order to minimise the background effects, 
an issue that will be discussed later in the paper. And it was made clear in Ref.(l4|] 
that the higher orders in the perturbative expansion in a influence considerably the 
differential decay rate, if the e + e~ -invariant-mass resolution becomes smaller than 
about 1% of the pion mass. This is so because, in that case, the first-order result 
already gives a larger-than-20% correction. The present experimental resolution 
for the e + e~ invariant mass being on the order of 0.5 MeV ~ 0.4% |15[| , an 
improved calculation of the radiative corrections for this decay is needed, in order 
to include higher orders in a. 

The following analysis is similar to the one usually applied to calculate the 
initial-state radiative corrections to meson production via electron-positron colli- 



sions [16]. The differential decay rates, denoted by P(x) below, where x is the 



invariant mass of either the pion or the e + e~ pair, are defined here to be equal to 
~dx 6 ^ • Assuming a relativistic Breit-Wigner shape for the pion resonance, 
the result for P{s) to zeroth order in the fine structure constant a is 

*(•) - *M) ( ._ m ff +Itmj . ( 2 > 

where yfs is the invariant mass of 7r°, m n = 134.9 MeV is its mass, and ~ 8 eV its 
natural width The pion natural width is so small, that any possible dependence of 
Po on s coming from the form factor associated to the decay ir° — ► e + e~ is neglected. 

Generally, if the experimentally measured quantity is the invariant mass "v/s 7 
of the e + e~ pair, the non-corrected and corrected differential decay rates, denoted 
by Pq(s) and P ra a( s ) respectively, can be both expressed as follows: 

f 2s ' 

^0,rad(s ) = J , dsP o(s)fo,v&d{s - S ) (3) 

The distribution fo,rad( s ~ s ) is subject to the requirement that 

J, ds/ ,rad(s - S ) = 1. (4) 



In a reference frame where the two final electrons have equal and opposite 3- 
momenta, the relation s — s = <im e u: holds, where uj is the total energy of the 
radiated photons, and m e is the mass of the electron. The upper bound (2s ) of the 
above integration is placed in order to avoid the pair creation of e + e~ pairs having 
invariant mass s and originating from the radiated photons, since they could be 



1 This is a different approach from the one in Ref. mM, where — ~7 e e , the quantity corre 



sponding to Po(s), is taken to be proportional to S(s — m 
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misidentified as being e + e pairs having the same invariant mass and originating 
from the pion. More on this upper bound will be discussed shortly. 

The differential decay rate Po(s') is associated to the distribution 

fo(s - s) = 6(s - s) (5) 

because, to zeroth order in a, s = s' . The effect of the radiation of soft photons 
is the "smearing" of this 6- function, even though the resulting differential decay 
rate still retains an integrable singularity at s = s . The convolution of the Breit- 
Wigner shape of Eq.(|2|) with the adequate "smearing" distribution / ra a(s — s ), the 
calculation of which is the main purpose of the present paper, gives us the radiatively 
corrected differential decay rate -P ra d(s ). 

The pion is assumed to be produced hadronically, as given, for instance, by 
the hadronic reaction K + — > and therefore only final-state QED radiative 

corrections are considered. Furthermore, it is assumed that all radiated photons, 
no matter what their energy is, escape detection. The calculation presented here 
is restricted to the soft-photon resummation technique, since this is known to give 
the bulk of the radiative corrections. The final result, applicable only to e + e~ 
invariant masses very close to the pion mass, is non-perturbative, since the leading 
logarithms of diagrams of all orders in a are summed, in a process well known as 



"exponentiation" |17 |. 



In this context, P rac j can be expressed as an infinite sum of individual dif- 
ferential decay rates P(ir° — > e + e~ + 7x7), each having n real soft photons radiated 
from the 7T°: 
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Prad(s') = E P ( 7r ° - e+e ~ + n 7)- (6) 
n=0 

Intuitively, this infinite sum has the meaning that in the soft-photon (classical) 
limit, where radiation loses its particle-like properties, two final states containing 
n and n + 1 photons respectively are undistinguishable, so a summation over all 
possible numbers of photons is required. In such an analysis, each of the quantities 
P(ir° — > e + e~ + nj) contain the virtual photon corrections needed to cancel the 
infrared divergencies. Moreover, the n soft photons are considered to be independant 
from each other. For this to be true, their emission must not influence considerably 
the motion of the electrons, so they must have an energy much smaller than the 
electron. Therefore, the relation s — s' <C 4nig must hold, and then the logarithms of 
the form ln(s/s — 1) dominate the QED perturbative expansion, and are summed 
to all orders in a. 

The appearance of these logarithms is closely connected to the cancellation 
of the infrared divergencies. In order to see this, it is assumed that the emission 
of soft real photons, of energy less than uq, is used in order to calcel the infrared 
divergencies associated to the virtual photons. Then the calculation involving soft 
real and virtual photons contains logarithms of the form In (^h^J ■ On the other 
hand, the hard-real-photon part of the calculation, involving photons of energy 
larger than uo, contains logarithms of the form In ( S Z? ). The addition of these 
two contributions cancels the dependance of the radiative corrections on u>o, leaving 
logarithms of the form In (s/s — 1) in the final result. 

The soft photons are assumed to be radiated only by the two real final 
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electrons. This is because the pion decay cannot proceed via the emission of three 
photons, because of C-conservation, so no photon can come from there. Moreover, 
the effect of photons emitted by the virtual electron, denoted by e* on the right side 
of the "77c*" triangle (see Fig.[lJ) , is neglected. Intuitively, this happens because 
the radiated soft-photons, which give the largest corrections, correspond to large 
spatial scales, and are therefore associated to the two final real electrons. The 
virtual electron, however, is far off-shell, since the typical momenta flowing in the 
triangle of FigJT] are on the order of m^, and it is associated to small spatial scales, 
so its contribution to these corrections is negligible. 

According to the previous discussion, only radiative corrections pertinent to 
the two final electrons are considered, and the result is, according to the analysis of 
Ref. 0@ 



fa-s'V 3 - 1 A 2 
/ rad ( S - s) = rJ—Al ad + -S* (7) 



with 



f3 = f3(s') = — (ln(s'/ml)-l). (8) 

This expression gives a reasonable estimate only in the relativistic limit s' 3> m 2 , a 
relation that obviously holds for s' near the pion mass. Moreover, this estimate is 
even more accurate when /?< 1. Here /^(m 2 ) ~ 0.047. 

1 2 

The quantities A J. correspond to radiative corrections that can be treated 
perturbatively. Among various other terms, they contain powers of /3, i.e. basically 
terms of the form (a In (s /m^)) n , that were left out of the exponentiation procedure. 
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The effects of vacuum polarization of the two virtual photons are also assumed to 
have been absorbed in A rad . To zeroth order in a, Aj ad = 1 and A^ ad = 0. In the 
following, the approximation 

A lad « 1 (9) 

and 

aL « o (io) 

is used, i.e. only the factor multiplying Aj ad , the non-perturbative part of the 
calculation, is considered, because it is expected to give the bulk of the radiative 
corrections. In view of the magnitude of /3, the error induced by this approximation 
should be on the order of 5%. 

The distribution / ra d( s — s ) contains the exponentiated logarithms, which 
give the leading QED correction to the decay profile for s « s . Since f3 < 1, 
/rad( s — s ) h as an integrable singularity at s = s' . It is worth noting that this form 
of /rad(s — s ) implies the assumption that the logarithms In (s/s — 1) dominate 
the QED perturbative expansion of the radiative corrections. This ceases to be 
true when, in the integral of Eq.(j3|), the integration variable s approaches the upper 
bound (2s') of the integration. Moreover, in that case, interference effects between 
this process and the main background process, which will be discussed shortly, 
become non-negligible. However, the present form of the "smearing" distribution 
/rad(s — s ), and the smallness of (5 in Eq.@, makes the integral of Eq.(||) completely 
dominated by the region of s ~ s . Therefore, effects coming from a large s — s 
difference are not expected to influence the final result considerably. 

In Fig. |, the quantity p 2 ( ^ ) -Po(g)/rad(g - s), evaluated at s = ml, is 
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Figure 2: The integrable singularity of the quantity p S , 1 ■2\ Po(s)fraxi(s — s), eval- 
uated at s = m^. 

plotted as a function of the pion invariant mass minus m n , i.e. y/s — m n . The 
singularity at s = s' shows clearly. The area under the curve that is partially 
shown in Fig.|2|, contained in the region < yfs — < , is equal to 

P rad (ml)/P (ml). 

It is not difficult to see from Eq.(^) and (g) that, for T w -C and (3 <C 1, 
P TaA (ml) » P (ml) (^-Y , (11) 



9 



with P evaluated at s = m^. The ratio appears naturally, as being the charac- 
teristic energy scale over which the pion differential decay rate varies considerably, 
divided by the largest energy scale of the process. The dependance of the right- 



hand side of Eq.(ll) on m e comes only through the parameter (3, which contains 
the term In (s'/mg), responsible for the collinear singularity in the limit m e — > 0. 
The fact that m e does not appear anywhere else in the final result is attributed to 
the factorisation of collinear singularities, a well-known property of QCD, which is 
also shared by QED [18]. For the given values of {3(m^), rV and m n , Eq.(ll) leads 



to the relation 

Prad(m^) « 0.46 P (ml). (12) 

Therefore, the peak value of the radiatively corrected differential decay rate 
is about 50% smaller than the one without radiative corrections 0. The magnitude 
of this effect shows that, for e + e~ invariant masses very close to m^, a calculation 
of these corrections using perturbation theory up to first order in a would not give 
an accurate estimate. 

The position of the maximum of -P r ad(<s') lies approximately at 

V7 Rj 77V - "^-IV, (13) 

o 

which corresponds to a shift of about 0.2 eV from m n , and is therefore experimen- 
tally unobservable. 

2 Similar examples can be found in connection with the production cross-section of heavy mesons 
via e + e~ colliding beams, where the QED radiative corrections are predicted to have a 50% effect 
Jl6| ]. In those cases, however, in contrast to the analysis performed in this paper, the characteristic 
energy scale over which the production cross-section varies considerably is determined not by the 
width of the produced resonance, but by the energy spread of the colliding beams, which is much 
larger than the width. 
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Figure 3: The lowest order and radiatively corrected differential decay rates Po(s ) 
(upper curve) and P ra ,d(s ) (lower curve) respectively, normalized over Po(m^), as 
functions of the electron pair invariant mass minus m n , i.e. (s') 1//2 — m n . 

The differential decay rates Po an d P ra ,d as functions of the invariant mass 
of the electron-positron pair minus m n , i.e. — m n , given by the numerical 

evaluation of the integral in Eq.(||), are shown in Fig. ||[ The effect of soft-photons 
radiated in the final state is apparent, not only in the maximum height of the 
radiatively corrected differential decay rate, but also in the radiative "tail" that 
appears for lower e + e~ invariant masses, as expected, since the photons are radiated 
in the final state. Indicatively, numerical results show that the radiatively corrected 
differential decay rate is roughly three times larger than the non-corrected one at 
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e + e~ invariant masses 150 eV below m n , and two times smaller at 150 eV above m w . 
In both cases, these differential decay rates are roughly three orders of magnitude 
smaller than the ones at s' ~ m%. This situation is to be contrasted to the case 
of initial-state radiative corrections, for processes like production of resonances via 
e + e~ collisions, where one sees a more pronounced radiative tail for invariant masses 
larger than the resonance mass. 

The experimental resolution of the invariant mass of the e + e~ pairs being 
on the order of 0.5 MeV, the structure of the decay curve shown in Fig.|| is not 
experimentally observable. By numerically integrating the function -P ra d( s ') over 
s , it is found that the number of e + e~ pairs detected in an energy bin of width 
0.5 MeV, centered around the pion mass, is approximately 80% of the number it 
would be without radiative corrections, i.e. we have a 20% correction. This happens 
because the radiative tail for s < partly compensates the reduced number of 
e + e~ pairs at s ^m^. A complete analysis of this process should take into account 
this binning for energy regions further away from m n . However, results obtained in 
this way should not be trusted for Vs too far away from (more than a few MeV), 
because of several reasons. First, because background effects become non-negligible. 
Moreover, in that case, one should add to this calculation the contribution of hard 
photons, which have an energy that is not small compared to the electron mass 
any more, and which become increasingly important for invariant masses further 
away from m n . Last, the dependance of the form factor, associated to the decay 
7T° —* e + e~, on s' cannot be neglected any more. 

It is furthermore expected that the effect of radiative corrections on the 
total number of e + e~ pairs, having all possible invariant masses, is dominated by 
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perturbative quantum effects, and that these are on the order of a few percent. 
In fact, a first-order in a prediction for this effect (Eq.(18) in Ref. |14|), which is 
a reasonable approximation, gives a 3% correction. It is therefore seen that the 
radiative tail of the differential decay rate, when considered throughout the whole 
range of allowable e + e _ invariant masses, compensates, to within a few percent, 
the reduction of the differential decay rate at s ~ m^. These radiative effects on 
the total number of e + e~ pairs are not going to be discussed further here, because 
they are completely overwhelmed by interference terms between this process and its 
background, and by the background itself, which is discussed here below, and they 
are therefore experimentally unobservable. 

To first order in the fine structure constant, the process considered here 
involves a final state of an electron-positron pair and at most one photon. A large 
background to this process is the same final state coming from the process 7r° — ► 77, 
and the subsequent decay of one of the final photons into an e + e _ pair flSfl - This 
process is of lower order in a, and it has a much larger decay rate than the one 
considered here. However, this background is easily distinguishable from the signal 
of interest, because its differential decay rate does not have a pronounced peak at 
e + e~ invariant masses near the pion mass. Moreover, the interference of this 
background with the process of interest in this paper is negligible in the region 
that is considered here, which is W m n (to within a few MeV) p0[ . 

To conclude, this paper has dealt with the shape of the ir° —* e + e~ differential 
decay rate, when radiation of soft-photons from the two final electrons is taken 
into account. In particular, the peak value of the radiatively corrected differential 
decay rate -P m d is smaller than the non-corrected one by about 50%. Moreover, 
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the radiative tail, emerging at invariant masses \fs smaller than m n , even though 
it has a considerably smaller height than the peak differential decay rate, it does 
not fall-off very rapidly with W, and, when considered over the whole range of 
allowable electron invariant masses, it almost compensates the reduction of the 
number of e + e~ pairs at W f» m n . This form of the differential decay rate should 
be taken into consideration when comparison between experimental and theoretical 
results for this decay is made, especially when stringent e + e~ invariant mass cuts 
are placed. A more involved analysis should include the effect of higher order terms 
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in the quantities A r ^ d . Finally, a similar analysis can be done for the decay of other 
mesons, like the ry°, to a lepton-antilepton pair. 
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